Polaron effects in electron channels on a helium film 
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Using the Feynman path-integral formalism we study the polaron effects in quantum wires above a 
liquid helium film. The electron interacts with two-dimensional (2D) surface phonons, i.e. ripplons, 
and is confined in one dimension (ID) by an harmonic potential. The obtained results are valid 
for arbitrary temperature (T), electron-phonon coupling strength (a), and lateral confinement (uJo). 
Analytical and numerical results are obtained for limiting cases of T, a, and ujq. We found the 
surprising result that reducing the electron motion from 2D to quasi- ID makes the self-trapping 
transition more continuous. 
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I. INTRODUCTION 

In recent years multi-channel electron systems have 
been realized over liquid helium by using one-dimensional 
(ID) metallic gate structures or by providing a surface 
distortion through a corrugated substrate in addition ta 
the holding electric field perpendicular to the surface.EJ 
A single channel was also obtained asing a polymer 
grooveB and strips on a circuit boardJj In these quasi- 
one-dimensional (QID) systems, the electron motion is 
quantized in the z direction by the effective holding po- 
tential due to an external electric field and image forces 
coming from the substrate, and in the y direction by a lat- 
eral confinement potential which can be modelled quite 
accurately, as in a wide class,of QID and QOD electron 
systems, by a parabolic well.a 

It is well known that the interaction of the electron 
with the quantized surface excitations (ripplons, in the 
case of a liquid surface) can be described as a po- 
laron problem with the amazing possibility of a phase- 
transition like behavior between localized and delocak 
ized states when varying the electron-ripplon coupling.El 
The ripplonic polaron state on the surface of helium re- 
sults into a decrease of the conductivity of the electron 
layer. More detailed information of the polaron state can 
be obtained by measuring the electron conductivity as a 
function of the microwave frequency.El Even though the 
formation of the localized polaron state has been claimed, 
there is no conclusive interpretation about the origin of 
the mobility dip iii some experiments and the question is 
still controversial.!] 

Quite recently ID confinement effects on surface elec- 



trons on bulk helium have been addressed by using the 
hydrodynamic model of the polaron which describes the 
energetics of the dimple formation on-the surface of liquid 
helium and its transport properties.! 

Even though surface excitations of the liquid helium 
film have in general a complicated dispersion relation 
coming from contributions of surface tension, gravity and 
film thickness, in the case of thin films the ripplon spec- 
trum has a well-defined acoustical character .0 In a previ- 
ous work,Q the 3D acoustical polaron was studied within 
the Feynman approach and it was found that, as a func- 
tion of the electron-phonon coupling constant a, the po- 
laron undergoes a self-trapping transition at a ~ \/fcoj 
where is a finite Debye cutoff in phonon space, which 
is continuous for fco < 18 and becomes discontinuous 
when fco > 18. For a 2D system, a similar self-trapping 
transition was found at a ~ 0.25 which is almost k^- 
independent. The transition is now continuous (discon- 
tinuous) for fco < 130 {ko > 130). E3 

In this paper, we will investigate how the polaron ef- 
fect is altered when we further confine the electron sys- 
tem to a channel, i.e. a QID system. In the case of 
LO-phonon interaction it is well known that the polaron 
effect (renormalization) increases with_.pQnfinement, i.e. 
with reduction of the dimensionality.Einlj Here, we will 
study the competition between confinement and polaron 
effects as resulting from the interaction with acoustical 
phonons (ripplons) . The confinement alters the available 
phase space of allowed fc- vectors for virtual phonon emis- 
sion and absorption and in this way the polaron effect. 

The present paper is organized as follows. In Sec. 
II, we present the Frohlich-type Hamiltonian, the path- 
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integral Feynman formalismJli'Q the trial Hamiltonian 
used and the variational principle to derive an upper 
bound to the exact free energy of the system studied. 
Explicit analytic results in limiting cases of temperature 
T, electron-phonon coupling strength a, and lateral con- 
finement ujQ are presented in Sec. III. The numerical 
results are shown in Sec. IV. In the last section we will 
present our conclusions. 



II. HAMILTONIAN MODEL AND FREE 
ENERGY 

The Hamiltonian which describes an electron interact- 
ing with 2D phonons and subject to a lateral confinement 
due to an harmonic potential is given by 



H 
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EK-(«K~+|) + ^/' (1) 



with the electron-phonon interaction given by 

Hi = Y. {^iH^'-' + y~>h-''-') . (2) 
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where r [p) is the 2D position (momentum) operator 
of the electron with mass m, Vj: is the electron-phonon 
strength and at (ar) are the creation (annihilation) oper- 

k 

ators of the phonons with wave vector fc, and frequency 
wj. The electron is laterally confined by an harmonic 
potential with frequency u)q. 

The Helmholtz free energy is given by 



F 



logZ, 



with 
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the partition function. The free energy, Eq. (H), and 
the partition function, Eq. (^), contain the contribu- 
tion of the free phonons, Fph- Here, we will shift the 
free energy F by this constant contribution Fph- Using 
the patb-iTit<?gral formalism, the trace of the partition 



functior 



can be written as 



dr I Dr{u)ex:p{S[r{u)]}5{r{P) - r)5(f(0) - r), 



(5) 

where / Dr denotes the integral over all possible electron 
paths, /3 = l/ksT is the inverse temperature, and u is 
related to the real time variable t as u = it. The action 
S[r\ in Eq. (^) is obtained after the exa^ elimination of 
the phonon coordinates and is given b}il3 

S^Se + Si, (6a) 



with 
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(6b) 



, (6c) 



where 

GUt) = ^[(1 + n{u;))e-''^\'\ + +7i(c.)e''-l*l] (6d) 

is the phonon Green's function, and n{Lu) = (e'^''" — 
the number of phonons with frequency u. 

The path- integral in Eq. (|^), can not be |-e.valuated 
exactly. We follow Feynman's polaron theoryllj and in- 
troduce a trial action Sq, such that: i) the path integral 
with this action can be done exactly, and (ii) it approxi- 
mates the original action S as close as possible. The free 
energy associated with this trial action is given by 

Zo = e-f^^" ^ j '^^ j ^Kw)exp{S'o[f(u)]} x 

5{m - r)5{m - r), (7) 

and we define the expectation value of any functional A[r\ 
with respect to 5*0 by 

{mo = ^ = / ^^»e^p{^o[^»]} 

A[r{u)]6{r{l3) - r)6{r{0) - r). (8) 

Using this trial action, the partition function of the 
real system, Eq. (|^), can be written as 



Z = Zo(exp [S - 5o))o. 



(9) 



If the actions S and Sq are real, as is the case for our 
problem, the convexity property of the function exp {x) 
can be applied from which we find the inequality 

(exp(5-5o))o>exp((5-5'o)o). (10) 
leading to a variational principle for the free energy 



(11) 



We propose a trial action derivable from the trial Her- 
mitian Hamiltonian 

^°=2;;^+2M + 2^("-^) +2;^+2"^" ^' 

(12) 

where X (px) is the ID position (momentum) of the 
fictitious particle of mass M, K is the strength of the 
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harmonic potential and u) is the oscillation frequency 
associated with the lateral confinement of the electron in 
the y direction. The first three terms in Eq. (^2|) corr|e=. 
spond to the 1-D Feynman polaron in the x dircction.tZI 
The last two are associated to the confined motion and 
the variables M, K and lo are variational parameters. 

Using the same procedure as for the original Hamilto- 
nian, Eq. (|I|), we eliminate the coordinates of the ficti- 
tious particle, and obtain the trial action 
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^ 1 „ /-/^ 

dt{x{tf + y(tf ) - -mJ \ dty(tf 



dt 



dsG^{t-s)ix{t)~x{s)f, (f3) 



with uj ^ y^K/M and = K[{M + m) /Mm]. 

Followmg the same approach used by Peeters and 
DevreesetS, from Eqs. ( |6a| ) and (|f3|), the Feynman vari- 
ational principle for the free energy of our system, Eq. 
fO), can be written as 



F<Ff-Fb.--(S~So)q, 



(14) 



where Fp is the free energy of the trial Hamiltonian, Eq. 
(Ill), given by 



Ff = -^log[L^-— ^] 
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2sinhi/3?; 2sinhi/3( 



xr7}], (15a) 



Fr is the free energy of the fictitious particle without the 
interaction term, given by 
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^«^~^^°^[2sinhi/3c.^ 



(15b) 



and 
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GUt-s){{x{t)-x{s)f). 
Defining 

{e''^<m-Pis))^^ ^A,ik,,t- s)Ay{ky,t - s) 
where 

A^{k^,t-s) = {e^f^dm-as))^^ 
with = x,y, we have that 



(15c) 



(16) 
(17) 



mt) - e(s))')o = -^A^{k^,t-s). (18) 
ok^ 



The three averages in Eq. (15c) can be calculated in 
an analogous way as done by Osakall3 and Peeters and 
DevreseeJla such that we obtain 



-(yW^>o-2^coth^, 



A^{k^, t — s) — e 



-ki:'^D((t~s) 



(19) 
(20) 



with 



Vt^ 



3 (l-e-^'l*! -4n(«)sinh^^), (21a) 

and Dy (t) is obtained from the above expression by tak- 
ing Lu and u — > w', giving 

Dy{t)^^{l-e-'^'^-4n{v)smh''^) . (21b) 



Substituting Eqs. (|l5a|)-(|21b|) in Eq. we finally 

obtain an upper bound F^ for the exact free energy of 
our system, given by 
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/ sinh ^I3uj 
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{u'' u;o') coth V - E l^fcl'tl + "K-)] (22) 
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(^yg-'^fc"g-fex D^{u)^~ky Dy{u)^ 

which is subject to the minimization conditions 

dF^ QF-" dF"" 



dv duj duj' 



0. 



(23) 



III. GROUND-STATE ENERGY AND LIMITING 
CASES 



The result obtained in Eq. (^2[) is general and can 
be used to calculate all thermodynamic quantities, like 
e.g. specific heat, entropy, internal energy, etc. In the 
zero-temperature limit, the free energy F'" reduces to the 
polaron ground-state energy Eq" ^ given by 



En''^ lim F" = -uj' - 

fS^oo 2 
-. JQ 



-{v-ujf- 

4:V 
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(24) 
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with 



and the polaron energy can be written as 



2w3 
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(25a) 



and Dy{t) is obtained from the above expression by tak- 
ing w — > and u — s- w', giving 

1 



(25b) 



In Eq. (|2J) the first term corresponds to the zero-point 
energy of the harmonic oscillator with frequency lj , the 
second one is due to the interaction between the lateral 
confinement and the fictitious oscillator with frequencies 
loq and w , respectively. The third term is due to the 
Feynman polaron with eigenfrequency v and the fictitious 
particle with frequency w (i.e. m — > oo). The last term 
is associated with the electron-phonon interaction. 

If the ^-confinement is much larger than the phonon 
effects, one can take the limit cj' w wq- Also, if in this 
limit the phonon cloud surrounding the electron is very 
small, we can assume v ^ w. Within these approxima- 
tions, the polaron ground-state energy E'o", Eq-(24), is 
given by 

t 1 

771 con J 



E 
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k P / (- 



(26) 



which equals the result from second-order perturbation 
theory. 

The first term in Eq'^°'^^ is the zero-point energy of 
the harmonic oscillator, the second exponential term un- 
der the integral is associate with the free motion along 
the X direction, and the last exponential under the in- 
tegral is due to the oscillator motion in the confinement 
direction containing all confinement levels. Expanding 
the last exponential term under the integral of Eq'^""^ , 
we can integrate out the u variable resulting in 

1 
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n=0 
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-ky /2loo 



(27) 



which allows a perturbation analysis in terms of dia- 
grams. 

We now consider the limiting cases of weak and strong 
coupling. For weak coupling, we have the following 
conditions: a 3> 0.5, v ~ w, and uj' 0. There- 
fore, expanding the exponentials in the D^.y terms, we 
obtain Dx{u) « u/2 and Dy{u) w u/2, resulting in 
f{k,u) ~ exp(— fc^u/2). In this case, the integral 

fko /• oc 

B^a I dkk^ I due-^''f{k,u) 
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(28) 



1 /cj. 



B. 



(29) 



In the absence of confinement [loq — 0), even when uj' is 
too small, Epoi = —B. For o^o > 1, the polaron energy 
increases fast. 

In the strong limit, a ^ 0.5, we have the condi- 
tions: V ^ w, and lu' v. Therefore, as before, ex- 
panding the exponentials in the Dx^y terms we obtain 
now Dx{u) ~ l/2v and Dy{u) « l/2a;', resulting in 
f{k,u) ~ exp(— fc^/2ti). In this case, the integral 



B = a I dkk^ J due'~''''f{k,u) 



~ -av (e-'=o/2- - 
and the energy is given as 



E,^ 



pol 



2 4u 



(30) 



(31) 



One can see that the energy decreases when we increase 
the lateral confinement. 



IV. NUMERICAL RESULTS 

In this section we present the numerical results for the 
zero-temperature limit of the free energy F^, that is, the 
confinement polaron ground-state energy Eq" , by min- 
imizing Eq. (p^ , and considering the electron interact- 
ing with longitudinal surface acoustic phonons (ripplons). 
For these excitations, the dispersion relation is given by 



s\k\ 



(32a) 



where s is the velocity of sound, and the Fourier trans- 
form of the electron-phonon interaction 



\Vk\ 



ak, (32b) 

describes both 2D phonons and ripplonsi0, where a 
is the dimensionless electron-phonon coupling constant 
which depends on the deformation potential. In our cal- 
culation the sum over the phonon wave vectors will 

be replaced by the integral A/ {2tt)'^ J dk, which is cut- 
off at fco: the Debye critical wave vector in phonon space 
which simulates the discreteness of the lattice, fco l/^i 
with a the lattice constant, and corresponds to the cap- 
illary constant kc ~ where d is the film thickness 
in case of liquid helium. Hereafter, we use dimension- 
less units and express the energy in units of ms^ and the 
length in units of h/ms. Substituting Eqs. ( |32aD -( p2h| ) 
in Eq. (|2^), and integrating out the angular coordinate 
of k we obtain 
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-a / dfcfc^ / due-'="/(fc,u), (33) 
Jo Jo 

with 

/o(^|i?."(u)-/?,"Hl). (34) 

The results for the energy as a function of a are shown 
in Fig. 1, for different values of the Debye cutoff ko and 
the lateral confinement frequency loq. The figure shows 
that the slope of the polaron energy varies at a ~ 0.51. 
This value of a is approximatelly twice the value obtained 
in Ref.lO, even when we take the limit ujq 0. This dif- 
ference can be understood as due to the fact that the 
translational invariance in the y direction is not taken 
into account in our model Hamiltonian, i.e, we are con- 
sidering a QID system. 

As can be seen, the lateral confinement increases the 
energy and this is more significant for small values of 
a. In order to analyze the effects of the lateral confine- 
ment on the self-trapping transition, Fig. 2 shows the first 
derivative of the ground-state energy as a function of a 
for two values of ko. We observe, in Fig. 2(a), that the 
lateral confinement, coq = 20, smoothens the curve com- 
pared with the one where uq — 0. Considering ko — 150, 
Fig. 2(b) shows the existence of a discontinuity in the first 
derivative corresponding to a first-order transition in the 
case where luq = 0. This is not observed for cuq = 20. 
These facts are best visualized in the second derivative 
of the ground-state energy as function of the coupling 
constant presented in Fig. 3. 

In Fig. 4, the phase diagram (/cq, i^o) for the 2D acous- 
tical polaron is depicted for a = 0.515 where the transi- 
tion occurs. For each loq, the value of fco determines the 
self-trapped transition and the curve bounds two regions. 
In the upper region, the transition is discontinuous, while 
in the below region it is continuous for this a-value. In 
the latter region it is possible to find a discontinous self- 
trapping transition when a > 0.515. 

In Fig. 5 we present the Feynman polaron mass M — 
(v/w)'^ as a function of a, for different values of the De- 
bye cutoff fco and the lateral confinement frequency luq. 
The figure shows clearly the transition from the weak- to 
strong-coupling regime. It also shows that for fco < 70 the 
transition is continuous. Again, it is worth to mention 
that, for ujQ = 0.0, the vaiue obtained here is different 
from the the 2D polaron.t2l Once more, this difference is 
due to the fact that the translational invariance in the y 
direction is not taken into account in our model Hamil- 
tonian. 

We also observe that the effects of the lateral confine- 
ment on the polaron mass are more significant in the re- 
gion of small a, i.e. the non self-trapped region, and the 
transition from low to large polaron mass is now smooth. 



In Fig. 6 we show the ground-state energy shift AE = 
E{ijjQ = Q) — E{uj{i) as a function of wq, for different values 
of fco and a fixed value of a. The lateral confinement in- 
creases the binding energy and we do not see any abrupt 
change in the energy values. This fact can be under- 
stood from Eq. (^7|). This result shows that, when cjo 
increases, the zero-point energy increases more than the 
terms corresponding to the free motion along the x di- 
rection, and the one associated to the oscillator motion 
in the confinement direction, that contain all the con- 
finement levels. Although the ground-state energy does 
not present an abrupt change with the confinement this 
fact is observed in the Feynman polaron mass. In Fig. 7 
we present the Feynman polaron mass M = (v/w)'^ as a 
function of ujq , for different values of fco and a fixed value 
of a. As can be seen, the polaron mass increases with the 
lateral confinement and can change dramatically for large 
values of fco where it exhibits a discontinuous behavior. 

V. CONCLUSIONS 

In this work we applied the Feynman path-integral 
method to study the effects of lateral confinement on 
the ground-state properties of the 2D acoustical polaron 
which can be formed in QID electron systems over a he- 
lium film. We evaluated the ground-state energy and the 
polaron mass as a function of the electron-phonon cou- 
pling a and we analyzed the effect of lateral confinement 
on these properties. We determined the phase diagram 
for polaron formation and found the intriguing result 
that the self-trapping transition becomes more contin- 
uous when the electron motion is reduced from 2D to a 
QID system. 
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FIGURES 

Fig. 1. Ground-state energy of the 2D acoustical polaron 
as a function of the electron-phonon coupling constant a 
for different values of fco, wq = 20 (dash-dotted lines) and 
Wo = (solid lines). 

Fig. 2. First derivative of the ground-state energy of 
the 2D acoustical polaron as a function of the electron- 
phonon coupling constant a for: a) = 50, b) ko = 150, 
in both cases for wq = O(solid lines) and wq = 20 (dash- 
dotted lines). 

Fig. 3. Second derivative of the ground-state energy of 
the 2D acoustical polaron as a function of the electron- 
phonon coupling constant a for: a) fco = 50, and b) = 
150, in both cases for uq = (solid lines) and lvq = 20 
(dash-dotted lines). 

Fig. 4. Phase diagram (fco, wq) for the self-trapping tran- 
sition of the acoustical polaron at a = 0.515. 

Fig. 5. Feynman polaron mass M = {v/w)'^ as a function 
of the electron-phonon coupling constant a for different 
values of fco and wq. 

Fig. 6. The ground-state energy shift of the polaron, 
Ai? = E{LOf) = 0) — E{ijj{)) as a function of the lateral 
confinement frequency ujq for a = 0.515, considering: 
fco = 50 the solid line, fco = 100 the dashed line, and 
fco = 100 the dash-dotted line. 

Fig. 7. Polaron mass as a function of the lateral con- 
finement frequency wofor the same parameters as in Fig. 
6. 
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